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ABSTRACT 

Molecular globules and pillars are spectacular features, found only in the interface region between a molecular cloud and an H n-region. Impacting 
far-ultraviolet (FUV) radiation creates photon-dominated regions (PDRs) on their surfaces that can be traced by typical cooling lines. With the 
GREAT receiver onboard SOFIA we mapped and spectrally resolved the [Cn] 158 /jm atomic fine-structure line and the highly excited l2 CO 
J=ll— »10 molecular line from three objects in Cygnus X (a pillar, a globule, and a strong IRAS source). We focus here on the globule and 
compare our data with existing Spitzer data and recent Herschel open-time PACS data. Extended [C n] emission and more compact CO-emission 
was found in the globule. We ascribe this emission mainly to an internal PDR, created by a possibly embedded star-cluster with at least one early 
B-star. However, external PDR emission caused by the excitation by the Cyg OB2 association cannot be fully excluded. The velocity-resolved [C n] 
emission traces the emission of PDR surfaces, possible rotation of the globule, and high-velocity outflowing gas. The globule shows a velocity 
shift of ~2 km s~ ' with respect to the expanding H n-region, which can be understood as the residual turbulence of the molecular cloud from which 
the globule arose. This scenario is compatible with recent numerical simulations that emphazise the effect of turbulence. It is remarkable that an 
isolated globule shows these strong dynamical features traced by the [Cll]-line, but it demands more observational studies to verify if there is 
indeed an embedded cluster of B-stars. 

Key words, interstellar medium: clouds - individual objects: Cygnus - molecules - kinematics and dynamics - Radio lines: ISM 
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Molecular pillars and globules are visually spectacular features 
that are produced when the ionizing radiation from OB-stars 
forms elongated structures in the neutral gas at the interface 
between Hn-regions and molecular clouds. Though known for 
a long time, these features were popularized by the Hubble 
Space Telescope (e.g. the famous 'Pillars of Creation' in M16). 
More recent Herschel observations from the HOBYS key pro- 
gram (Motte et al. I2010I I have revealed many nice examples 
of pillars and globules in Cygnus-X, our SOFIA (Stratospheric 
Observatory For Infrared Astronomy) target. An important ques- 
tion is under which conditions which type of stars can form in- 
side pillars and globules. So far, low-mass stars have been found 
(Sugitani et al. l2002j i but to date no clear observational signature 
of more massive stars exists. In this Jetterwe show observations 
that point toward the existence of a cluster of B-stars associated 
with a globule. 

Pillars have a column-like shape and a physical connec- 
tion to the gas reservoir of the molecular cloud while glob- 
ules are isolated and have a head-tail structure pointing toward 
the illuminating source (see Fig. [TJ. Both are only found in re- 
gions of OB-stars where a direct impact of U V-radiation creates 
photon-dominated regions (PDRs) on the surface of a molec- 
ular cloud/clump, often visible as a bright rim at the edge of 
the cloud. These PDRs are best traced in their cooling lines, i.e. 
atomic fine-structure lines of ionized carbon ([Cn] at 158 yum) 
and atomic oxygen ([Oi] at 63 and 145 yum), as well as high-J 
CO rotational lines. 
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Fig. 1. Three-color (3.6, 8, and 24-fim) Spitzer image of the 
southern Cygnus-X region, obtained from the Cygnus-X legacy 
survey (Hora et al. 12009, http://www.cfa.harvard.edu/cygnusX). 
The Cyg OB2 association is located at a projected distance of 
~1.2° north («40 pc, assuming a distance of 1.4—1.6 kpc, Rygl 
et al. 12012b . The sources observed with SOFIA (pillar, globule, 
IRAS20306+4005) are as indicated. 



The basic explanation for the formation of globules (e.g. 
Lefloch & Lazareff 1994) involves U V-radiation impacting on 



N. Schneider: Globules and Pillars seen with SOFIA 



a pre-existing clumpy molecular cloud and photoevaporating 
the lower density gas, leaving only the densest cores, which 
may collapse to form stars ('radiative driven implosion' sce- 
nario, Bertoldi 1989)- Recent (magneto)-hydrodynamic simula- 
tions including radiation (Gritschneder et al. 120091 Tremblin et 
al. l2012all2012bt begin to successfully model pillar/globule for- 
mation from a different perspective, emphasizing the importance 
of turbulence. 

To understand pillar and globule formation, in particular un- 
der which conditions stars can form within them, we initiated a 
large observational and numerical study. This program collects 
observations from a Herschel open time project], ground-based 
molecular line Mopra observations^], and the SOFIA observa- 
tions presented here. We use the HERACLES code (Gonzales 
et al. 120071 1 and th e radiative transfer PDR-code KOSMA-tau 
(Rollig et al. l20061 > for modeling. 

The initial objective of the [C n] 158 /mi and CO J= 1 1 -> 10 
SOFIA observations was to obtain spectrally resolved data to 
study the effect of external UV-radiation on pillars and globules, 
and to interpret these in light of numerical simulations. We chose 
to focus on the southern part of the Cygnus-X region (Reipurth 
& Schneider 2008 ) where the very massive and rich Cyg OB2 
association (with an estimated 2600 OB-stars, Knodlseder 2000) 
illuminates the molecular cloud (see Fig. [TJ. 



2. Observations 

The [Cn] atomic fine-structure line at 1.90 THz (158 /mi) and 
the CO J=l 1 — > 10 molecular rotation line at 1 .267 THz were ob- 
served with the Pi-heterodyne receiver GREAlQ on SOFIA on 
a flight 2011, November 10 from Palmdale, California. On-the- 
fly maps with a scanning speed of 6"/s and a size of 4'x3' and 
3.5'xl' were made of the globule and the pillar, and the IRAS 
source, respectively. The globule was taken in two coverages, 
while the pillar and the IRAS source were observed in only one 
coverage. The pillar does not have a sufficient signal-to-noise 
(S/N) ratio, and the IRAS source is not fully covered, there- 
fore we did not attempt a more quantitative analysis of these 
sources. We only use a spectrum of each (Fig. [6} here, because 
they assist our argument for globule formation. For blank sky- 
subtraction, an OFF-positions was used, based on molecular line 
and IR-surveys. A main-beam efficiency of 0.51 (L2 channel 
with [Cn] ) and 0.54 (LI channel with CO) was applied to the 
data (Heyminck et al. 2012) and a linear baseline was removed. 
The observed rms level is 3.3 K (on a main beam temperature 
scale in a channel of 0.23 km s~'), corresponding to a S/N ratio 
of « 6 for the peak position. The absolute calibration uncertainty 
is estimated to be around 20 % (Heyminck et al. 2012). 

In addition, [C n] 158 /mi, [Oi] 145 //m and 63 /mi data from 
the aforementioned Herschel open time project were obtained 
for the globule on 2010, December 14, using the PACS spec- 
trometer (Poglitsch et al. 2010). These observations consist of 
a single pointing (field of view 47"x47") in chopping/nodding 
range scan mode, with a chopper throw of 6' and Nyquist spec- 
tral sampling. The data were reduced using the standard reduc- 
tion pipeline in HIPE 7.0. The final maps, from which the line 
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1 'Pillars of creation: physical origin and connection to star forma- 
tion', PI: N. Schneider (PACS, SPIRE, and HIFI spectroscopy) 

2 Mopra Pis: P. Tremblin, V. Minier (CO, HCN, HCO + lines) 

3 The German REceiver for Astronomy at Terahertz frequencies. 
GREAT is a development by the MPI fur Radioastronomie and 
the KOSMA/Universitat zu Koln, in cooperation with the MPI fur 
Sonnensystemforschung and the DLR Institut fiir Planetenforschung. 
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Fig. 2. Top: Map of [C n] emission (color-coded, integrated be- 
tween 3 and 12 km s _1 in K kms 1 on a T m \, scale). The Op- 
position is RA(2000)=20 /, 33'"53 i Dec(2000)=40°8'45". The 
contours (1 to 3.5 K kms _I by 0.5 K kms -1 ) indicate the very 
confined CO J=ll— >10 emission. Bottom: Overlay of [Cn] 
contours (3 to 1 1 K kms~' by 2 K kms -1 in a T m b scale) in black 
on the Spitzer 5.8 /mi continuum emission. The central cluster 
and a smaller sub-cluster in the northwest in the head of the glob- 
ule (in pink) are clearly visible. The dashed rectangle indicates 
the area mapped in [C n]. Labels A, B, and the cross indicate the 
positions of the [C n] spectra shown in Fig. 4. 



fluxes are measured, were made with the PACSman software 
(Lebouteiller et al. 2012). These PACS maps (11.5" at 158 /mi) 
were then convolved to the approximate beamsize of SOFIA 
at 158 /mi (16") and the intensity ratio of individual positions 
was determined. The SOFIA [Cn] data are typically a factor 
1.6 stronger than the PACS data. Note, however, that PACS cal- 
ibration was only derived from on-source data (the off-source 
position was contaminated by emission), and the calibration for 
SOFIA is still being tested. Despite this, the peaks of emission 
coincide very well, within 5" (the PACS pointing accuracy is 
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Fig. 3. Overlay of PDR-tracer lines on the Spitzer 5.8 //in con- 
tinuum emission. White dashed contours: [Cn] 158 /im levels 
3 to 11 K kms -1 by 2 K kms" 1 . Green contours: CO J=ll->10 
levels 1 to 3.5 Kkms" 1 by 0.5 Kkms" 1 . Red contours: [Oi] 145 
/zm0.7e-16tol.5e-16by0.2e-16W/m 2 . 



~2"), indicating that the pointing of SOFIA, using planets and 
optical pointing with a guide-camera, has been very reliable. 
Table 1 summarizes the flux values of each target from both sets 
of observations. 



Center ~7.7 km/s 
A ~9.3 km/s J~ 



B ~6.1 km/s 




Fig. 4. Individual [Cn] spectra from the globule (angular res- 
olution 16"), with positions A, B, and center (+) as indicated 
in Fig. 2 (top). The velocities quoted in the panel are the line- 
center velocities obtained by a Gaussian fit. The 'center' spec- 
trum has the broadest linewidth (6.7 km s -1 ) compared to 4.0 
and 4.9 km s~' for positions 'A' and 'B', 



3. Results and analysis 

The [C u] map of the globule and an overlay on Spitzer 5.8pm 
emission is given in Figs. [2] and [3j whilst the [C n] spectra at 
three selected locations are found in Fig. |4] The [C n] emission 
covers the head of the globule with an increase in intensity to- 
ward the center. The CO J=l 1— »10 emission is very confined in 
the southwest as is also seen in Fig. [3j showing a zoom of the 
globule head and where we additionally plot contour levels of 
Herschel PACS [Oi] 145 pm emission. All PDR lines tracing in 
particular warm gas are confined in the same position. We over- 
lay [C n] and Spitzer on a near-IR 2MASS image in Fig. 7 where 



both images reveal a number of compact sources that are located 
at the peak of [C n] and CO emission. This is a cluster with at 
least one early-type B-star (see appendix A for details). To what 
extent this cluster is embedded has a strong implication on the 
interpretation of the [C n] emission distribution. 
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Fig. 5. Close-up of the [C n] map of the globule with the peak 
line velocity in color, determined from a Gaussian fit to each 
spectrum. From each Gaussian fit, the excess emission in the 
blue (~4 to 7 km s" 1 ) and redshifted (~9.5 to 11 km s -1 ) wing 
was determined and is plotted as blue and red contours (for better 
visibility the maps are resampled on a finer grid). The black tri- 
angles indicate the cluster memebers within the H n-region (note 
that this is where the hot gas, visible in CO J=l 1— >10 emission 
is found). The gray triangles mark the three sources from the 
secondary cluster. 



Though it is clear, considering the orientation and shape of 
the globule, that the impact of the Cyg OB2 association («40 
pc projected distance) created the globule, it is less obvious that 
the observed [C n] emission is only caused by externally illu- 
minated PDRs. In this case, a limb-brightened, arc-like emission 
distribution with a decrease of intensity toward the center should 
be seen. In contrast, [C n] [O i] and CO peak where the brightest 
sources are found, arguing for internal heating due to the clus- 
ter. However, because we do see 2MASS sources with IR excess 
projected just outside the globule, it is possible that the globule 
was larger in the past, formed an extended cluster, and B-stars 
from this cluster are now outside the shrinking globule, illumi- 
nating it externally. This may explain why the [C n] emission is 
quite extended and not as confined as the CO and [O i] emission. 
On the other hand, beam dilution and exitation by UV-photons 
leaking out to the surface with clumps having a small area-filling 
factor might also serve as an explanation. 

Studying the dynamics of the gas with the [Cn] spectra 
(Fig. |H, however, supports the idea that at least most of the 
B-cluster is embedded. The 'center', located at the position of 
the cluster, represents the bulk emission of the globule around 
7.7kms~ 1 with a broad linewidth, suggesting that material is 
moving toward and away from us along the line-of-sight (LOS). 
The spectra of positions A and B show a symmetric - but oppo- 
site - velocity shift of 1.6-2 kms -1 along the LOS. Both spectra 
show high-velocity wing emission, i.e. a prominent red wing for 
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the blueshifted 'B' spectrum and the reverse situation for the 'A' 
position. 

The velocity field of the central region of the globule is 
shown in Fig. 5. There is a sharp border between the two veloc- 
ity ranges (~6 and 9 kms -1 ), and the star-cluster is spatially lo- 
cated between these two velocity components. The velocity field 
shows a possible pattern of rotation around a slightly inclined 
northsouth axis of the globule, similar to what was already seen 
in CO-observations in pillars (Gahm et al. 120061 Hily-Blant et al. 
12005b . In addition, the globule is somewhat flattened orthogonal 
to the rotation axis, which could be interpreted as an additional 
argument for rotation. The minimum mass of the globule needed 
to support rotation is estimated from M = (v 2 r)/G) to 93 M 
with a velocity v of 2 km/s and radius r 0. 1 pc, which agrees 
with mass estimates from dust observations (see appendix A). 
Overlaid on the velocity field are blue and red contours of line- 
integrated [C n] emission from the two high-velocity wings, re- 
sulting in a typical 'outflow' pattern with two lobes. We attribute 
the high- velocity emission to warm gas that is accelerated by the 
expanding H n-region into an inhomogeneous environment. This 
scenario is similar to the one found in SOFIA [C n] observations 
of S106 (Simon et al. l20T2l 

The most probable explanation of the extended [C n] emis- 
sion and the confined CO and [O i] emission is therefore a com- 
bination of internal heating by the cluster, dumpiness in the 
globule, and some contribution from the FUV flux from Cyg 
OB2. To quantify this scenario, we determined an FUV intensity 
of x ~390 at the position of the globule head using the IRAS 
60 and 100 //m fluxes (Nakagawa et al. 119981) . External illumi- 
nation from Cyg OB2 can contribute to the flux, we estimated a 
range between^ a 1 80-370, depending on the number of O-stars 
(see appendix B for details) but is unlikely. In contrast, internal 
illumination by at least one B-star produces a flux of x ~360 
to 11000 (depending on the spectral type, see appendix A) in a 
distance of 0. 1 pc from the star, which would be consistent with 
our observations. 



Flux 
[10- 16 Wnr 2 ] 



Flux 

[10~ 5 erg s _I cirr 2 sr~'] 



4.34 


7.54 


0.7 (1.6) 


1.21 (2.78) 


2.71 


4.71 


1.24 (1.34) 


2.15 (2.32) 


2.34(2.51) 


4.07 (4.36) 



SOFIA [C II] 158 //m 
SOFIA 12 CO 11->10 
PACS[Cll] 158 yum 
PACS[Ol] 145 pm 
PACS[Oi] 63 pm 

Table 1. Flux values from SOFIA and Herschel PACS observa- 
tions of the globule. For comparison purposes, all maps have the 
same angular resolution of 16". The fluxes were determined for 
the position -40",0 (peak of [C n]) though for some lines, the 
value of the peak position is offset from this position (see Fig. 
3and is given in parenthesis. 



~2 kms _I with respect to the other two sources. The pillar and 
the IRAS source are still attached to the molecular cloud (see 
Fig. [TJ but directly face the expanding H n -region. We therefore 
define their velocity of ^6 km s -1 as that of the expanding H n - 
regioiQ. The velocity difference seen for the globule can then 
be understood from simulations (Tremblin et al., 12012b) as the 
residual turbulence of the molecular cloud. When the ionized 
gas pressure dominates the ram pressure of the initial turbulence 
in the molecular cloud, the dynamic is imposed by the expan- 
sion of the H ii -region, leading to a global expansion, with the 
formation of clumps and pillars in the dense border-shell of the 
primordial molecular cloud. When the ram pressure of the tur- 
bulence exceeds the ionized gas pressure, the cold gas has suffi- 
cient kinetic energy to enter the H n -region with motions that are 
perpendicular to the ionization propagation. Then more globules 
are formed that show a typical velocity shift with respect to the 
cloud of a few km s -1 . The shift can be regarded as the remnant 
signature of turbulence inside the molecular cloud. This picture 
is quite different from a radiative-driven implosion scenario of 
isolated clumps, which does not predict any particular velocity 
difference. 



IRAS ~4/6 km/s 
Globule ~8 km/s 
Pillar ~6 km/s 




Fig. 6. Representative [C n] spectra from the three observed re- 
gions, the pillar, globule (center position with bulk emission) 
and IRAS20306+4005. The double-peaked spectrum from the 
IRAS source is due to two seperate components at around 4 and 
6 km s -1 . The [C n] spectrum of the globule (red) shows a shift 
of around 2 km s -1 with respect to the pillar (blue) and IRAS 
source (black). 



4. Formation of the globule 

Comparing the [C n] spectra from the globule, pillar, and IRAS 
source (Fig.|6]l shows that the globule has a velocity difference of 
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Appendix A: The IR cluster 

Figure 7 shows a near-IR image from 2MASS|j with contours of 
Spitzer and [C n] overlaid. There is a higher stellar density at the 
location of the [Cn] emission, coinciding well in size, shape, 
and position. There is substantial crowding in the 2MASS im- 
ages and most of the fainter stars are typically red. On the red 
optical plates from the Digital Sky Survey only a handful of stars 
are seen. The extended emission seems brighter in the K-band. 
Judging from a preliminary J-H vs. H-K diagram, limited by 
source-confusion and sensitivity, the central cluster has a higher 
fraction of near-IR excess sources compared to the surrounding 
field. A preliminary color-color diagram also shows that sources 
with very different extinction values in the central cluster are 
found, though the median value is about 10 mag of visual extinc- 
tion of the 2MASS sources measured. Many of the faint sources 
detected by 2MASS have no reliable flux measurement, how- 
ever, and to make a detailed study of this cluster we need higher 
spatial resolution and deeper JHK images. 

Based on these 2MASS images, both Kronberger et al. 
(2006) and Kumar et al. ( 12006b classified the crowding of stars as 
' an associated, partly embedded cluster' . The latter give 9 cluster 
members within an effective radius of 0.53 pc, and a stellar mass 
of 17 M . Optical spectroscopy (Cohen et al. 11989b identified a 
B3 III star as one cluster member, but the spectral type classi- 
fication is very uncertain because it comes from a rough fit to 
Balmer line depths, and the luminosity class III is derived from 
the required flux and not from a proper spectral classification. 

Sridharan et al. (12002) derived a flux of 25 mJy at 3.6 cm and 
1 .4 Jy at 1.1 mm for this star, which is resolved at both wave- 
lengths, and they interpreted it as a 90 M clump comprising a 
high-mass protostellar object and an H n-region. Independently, 
the spectral type of the ionizing source can be estimated from 
the relation between the radio continuum flux density and the 
Lyman photon flux (Martin-Hernandez et al. 12005b . A flux of 
~13.7 mJy at 1.4 GHz has been measured (Setia Gunawan et al. 
12003b . which is consistent - as is the flux information at 3.6 cm 
- with an ionizing photon flux of F=1.8xl0 45 ph s , i.e. a Bl 
ZAMS star according to Panagia (11973b . 

Low-angular resolution millimeter continuum observations 
(BLAST, Roy et al. lDTTl SCUBA, Williams et al. 2005) indicate 
a dust temperature of 38 K and a lower limit of the total mass of 
about 40 M . Williams et al. determined a B2 spectral type from 
a luminosity of 4000 L and dust modeling. Using data from 
recent Herschel imaging within the HOBYS program (Motte et 
al. l2010b . we obtained a total mass of 80 M for the globule (size 
scale ^0.5 pc). 

Though the spectral classification of the internal B-star is 
fairly uncertain (see above), we can estimate a flux considering 
different spectral types. Far-ultraviolet fluxes for a number of 
early B-type stars, which include a Vir(BlIV), ceEri (B3V), and 
/? Cen (B1III), were measured by Holberg et al. (1982) using the 
UV spectrometers onboard the Voyager interplanetary probes. 
The spectral types of those three stars bracket the one estimated 
for the star at the center of the cluster, and their FUV luminosi- 
ties can thus be taken as upper or lower limits to the U V radiation 
internally injected in the globule. We used the distances to these 
stars as determined by Hipparcos (van Leeu wen 12007b to obtain 
their luminosities in the 912 A - 1225 A range. We obtained lu- 



minosities of l.lxlO 48 ph/s (a Vir), 9.7xl0 46 ph/s (a Eri), and 
3. Ox 10 48 ph/s (fi Cen). The corresponding volume-averaged val- 
ues of x, adopting a typical radius of 0.2 pc for the globule, are 
;tf=4,000 (a Vir), ;f=360 (a Eri), and ^=11,000 (J3 Cen). These 
values are higher, or at least comparable for the lower limit set 
by the B3V star a Eri, than the peak value of x derived from 
the IRAS 60 and 100 micron fluxes, thus supporting an internal 
origin for the excitation of the PDR. 

Appendix B: Estimation of UV-flux from Cyg OB2 

We here consider only O-stars because they dominate the UV 
emission over B-stars, with 06 stars dominating the overall 
emission output from the cluster. The number of O-stars in Cyg 
OB2 is very uncertain and estimates range between ~45 and 
~100 (see Reipurth & Schneider [2008 1 for an overview). In our 
estimates, we assume that Go = 10 6 at a distance of 0.1 pc from 
an O-star. The projected distance from the center of the OB as- 
sociation to the globule is 40 pc. Hence, 100 O-stars produce a 
flux of Go =* 10 8 /400 2 = 625 at the position of the globule, cor- 
responding to x — 370. We consider this a strict upper limit. For 
50 O-stars, this reduces to Go — 313 or x — 183. 

The flux from the IRAS 60 and 100 pxa. flux ratio is G = 660 
o^X =a 390 (obtained with a method described in Nakagawa et al. 
1998 ) in a large beam at the position of the globule, which agrees 
well with the values estimated above. We believe, however, that 
the contribution to the FUV flux from the Cyg OB2 association 
is significantly lower than the upper limit derived above due to 
the reasons detailed below. 

Because the volume around the OB association is not empty, 
there is still extinction on the way to the globule, which then does 
not receive the full UV flux from the entire O-star population. 
Unfortunately, it is not possible to give a precise quantification 
of the effect of extinction, but a reduction of the estimate of the 
flux by a factor of a few is possible. 

A basic problem in estimating the number of O stars is the as- 
sumption of a single age for the whole association. Knodlseder 
(2000) used 3 Myr and thus assumed that all stars are on the 
main sequence, with the brightest stars, therefore, being hottest. 
Hanson ( 120031) already showed that some of the stars are bright 
because they are (B-type) giants and supergiants, not because 
they are hot. There are currently around 70 spectroscopically 
confirmed O stars in the association, but more than half of them 
are 08-09. 

Although it is already difficult to specify a distance between 
the O-stars and the globule, the distance of 40 pc to the center 
of the association is a projected distance. Taking into account 
the depth along the line of sight can only increase the actual 
distance, and because the external flux decreases quadratic ally 
with the distance to the cluster, this in addition may significantly 
decrease the FUV flux received by the globule. 



5 The Two Micron All Sky Survey (2MASS) is a joint project 
of the University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by the 
National Aeronautics and Space Administration and the National 
Science Foundation. 



N. Schneider: Globules and Pillars seen with SOFIA, Online Material p 2 




Fig. 7. Near-IR image of the 6'x6' field centered on the cluster 
in the globule (RGB coded from 2MASS Ks, H and J) with con- 
tours (28, 30, 35, 50, 125, 200, and 600 MJy/sr) of the Spitzer5.8 
fim image (top) and the [C n] (bottom) contours (3-11 K km/s in 
steps of 2 K km/s) overlaid on it. North is up, east is left. 



